A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

02 (08.19.5 RX J1713.7-3946; 13.07.2) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Evidence for TeV gamma-ray emission from the shell type 
SNR RX J1713.7— 3946 



o 
o 
o 

(N 

c 

in 



> 

^- 

o 

o 
o 
o 

Of 

6 

CO 



X 



H. Muraishi 1 , T. Tanimori 2 , S. Yanagita 1 , T. Yoshida 1 , M. Moriya 2 , T. Kifune 3 , S. A. Dazeley 4 , 
P. G. Edwards 5 , S. Gunji 6 , S. Hara 2 , T. Hara 7 , A. Kawachi 3 , H. Kubo 2 , Y. Matsubara 8 , Y. Mizumoto 9 , 
M. Mori 3 , Y. Muraki 8 , T. Naito 7 , K. Nishijima 10 , J. R. Patterson 4 , G. P. Rowell 3 4 , T. Sako 811 , 
K. Sakurazawa 2 , R. Susukita 12 , T. Tamura 13 , and T. Yoshikoshi 3 

Faculty of Science, Ibaraki University, Mito, Ibaraki 310-8521, Japan 

Department of Physics, Tokyo Institute of Technology, Meguro, Tokyo 152-8551, Japan 

3 Institute for Cosmic Ray Research, University of Tokyo, Tanashi, Tokyo 188-8502, Japan 

4 Department of Physics and Mathematical Physics, University of Adelaide, South Australia 5005, Australia 

5 Institute of Space and Astronautical Science, Sagamihara, Kanagawa 229-8510, Japan 

6 Department of Physics, Yamagata University, Yamagata 990-8560, Japan 

7 Faculty of Management Information, Yamanashi Gakuin University, Kofu, Yamanashi 400-8575, Japan 

8 STE Laboratory, Nagoya University, Nagoya, Aichi 464-8601, Japan 

9 National Astronomical Observatory, Tokyo 181-8588, Japan 

10 Department of Physics, Tokai University, Hiratsuka, Kanagawa 259-1292, Japan 

11 LPNHE, Ecole Polytechnique. Palaiseau CEDEX 91128, France 

12 Computational Science Laboratory, Institute of Physical and Chemical Research, Wako, Saitama 351-0198, Japan 

13 Faculty of Engineering, Kanagawa University, Yokohama, Kanagawa 221-8686, Japan 



Received 1 December 1999 / Accepted 23 December 1999 



Abstract. We report the results of TeV gamma-ray 
observations of the shell type SNR RX J1713.7— 3946 
(G347.3— 0.5). The discovery of strong non-thermal X- 
ray emission from the northwest part of the remnant 
strongly suggests the existence of electrons with ener- 
gies up to ~ 100 TeV in the remnant, making the SNR 
a good candidate TeV gamma-ray source. We observed 
RXJ1713. 7-3946 from May to August 1998 with the 
CANGAROO 3.8m atmospheric imaging Cerenkov tele- 
scope and obtained evidence for TeV gamma-ray emission 
from the NW rim of the remnant with the significance of 
~ 5.6cr. The observed TeV gamma-ray flux from the NW 
rim region was estimated to be (5.3 ± 0.9[statistical] ± 
1.6 [systematic]) x 10~ 12 photons cm -2 s _1 at energies > 
1.8 ± 0.9 TeV. The data indicate that the emitting region 
is much broader than the point spread function of our tele- 
scope. The extent of the emission is consistent with that 
of hard X-rays observed by ASCA. This TeV gamma-ray 
emission can be attributed to the Inverse Compton scat- 
tering of the Cosmic Microwave Background Radiation by 
shock accelerated ultra-relativistic electrons. Under this 
assumption, a rather low magnetic field of ~ 11/xG is de- 
duced for the remnant from our observation. 



1. Introduction 

Supernova remnants (SNRs) are currently believed to be 
a major source of galactic cosmic rays (GCRs) from the 
arguments of energetics, shock acceleration mechanisms 
(Blandford & Eichler |1987| , Jones & Ellison |1991| ), and 
the elemental abundances in the source of GCR (Yanagita 
et al. |i"990| , Yanagita & Nomoto |i"999| ). EGRET observa- 
tions suggest that the acceleration sites of GCRs at GeV 
energies are SNRs (Esposito et al. 1996 ). However, direct 
evidence for the SN origi n of G CRs at Te V ener gies is 
scarce (e. g. Ko yama et al. 1995| , Allen et al. 1997 , Buck- 



ley et al. |199q ). Arguably the best evidence for the ex- 



istence of relativistic electrons with energies around 100 
TeV is the CANGAROO observation of TeV gamma-rays 
from the northeast rim of SN1006, which coincides with 
the region of maximum flux in the 2-10 keV band of the 
ASCA data (Tanimori et al. |l998b|). This TeV gamma- 
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ray emission was explained as arising from 2.7 K Cosmic 
Microwave Background Radiation (CMBR) photons be- 
ing Inverse Compton (IC) up-scattered by electrons with 
energies up to ~ 100 TeV and allowed, together with the 
observation of non-thermal radio and X-ray emission, the 
estimation of the physical parameters of the remnant, such 
as the magnetic field strength (Pohl [1996 , Mastichiadis 
1996|, Mastichiadis & de Jager 1996[ Yoshida & Yanagita 



1997, Naito et al. 1999) 
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The shell type SNR RX J1713.7— 3946 was discov- 
ered in the ROSAT All-Sky Survey (Pfeffermann & 
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Aschenbach 1996). The remnant has a slightly elliptical 
shape with a maximum extent of ~ 70'. The 0.1-2.4 keV 
X-ray flux from the whole remnant is ~ 4.4 x 10~ 10 erg 
cm~ 2 s _1 ranking it among the brightest galactic super- 
nova remnants. Subsequent observations of this remnant 
by the ASCA Galactic Plane Survey revealed strong non- 
thermal hard X-ray emission from the northwest (NW) 
rim of the remnant that is three times brighter than that 



from SN1006 (Koyama et al. 1997). The non-thermal emis- 



sion from the NW rim dominates the X-ray emission from 
RX J1713.7— 3946, and the SNR X-ray emission as a whole 
is domin ated by non-thermal emission (Slane et al. 1999 , 
Tomida 1999 ). It is notable that the observed emission 
region of hard X-ra ys ext ends over an area ~ 0°.4 in di- 
ameter. Slane et al. ( 1999j) carried out 843 MHz radio ob- 
servations using the Molonglo Observatory Synthesis Tele- 
scope, and discovered faint emission which extends along 
most of the SNR perimeter, with the most distinct emis- 
sion from the region bright in X-rays. Slane et al. (1999) 
suggest the distance to RX J1713.7— 3946 is about 6 kpc 
based upon the observation of CO emission from molecular 
clouds which are likely to be associated with the remnant. 
The dominance of non-thermal emission from the shell 



is reminiscent of SN1006. Koyama et al. (1997) proposed 
from the global similarity of the new remnant to SN1006 
in its shell type morphology, the non-thermal nature of the 
X-ray emission, and apparent lack of central engine like a 
pulsar, that RX J1713.7— 3946 is the second example, af- 
ter SN1006, of synchrotron X-ray radiation from a shell 
type SNR. These findings from X-ray observations would 
suggest that TeV gamma-ray emission could be expected, 
as observed in SN1006, from regions in the remnant ex- 
tended over an area larger than the point spread function 
of a typical imaging telescope (<~ 0°.2). 

Both SN1006 and RX J1713.7— 3946 show notably 
lower radio flux densities and relatively lower matter den- 
sities in their ambient space when compared with those 
for the other shell type SNRs ( Green 1998 ) for which the 
Whipple group (Buc kley et al. |1998|) and CANGAROO 
group (Rowell et al. 1999 ) have reported upper limits to 
the TeV gamma-ray emission. These characteristics might 
be related to the reason why TeV gamma-rays have been 
detected only for SN1006 and not from other shell type 
SNRs: the lower radio flux may indicate a weaker mag- 
netic field which may result in a higher electron energies 
due to reduced synchrotron losses. In addition, the lower 
matter density would suppress the production of ir° de- 
cay gamma-rays. An observation of TeV gamma-rays from 
RX J1713.7— 3946 would provide not only further direct 
evidence for the existence of very high energy electrons 
accelerated in the remnant but also other important infor- 
mation on some physical parameters such as the strength 
of the magnetic field which are relevant to the particle 
acceleration phenomena occurring in the remnant, and 
would also help clarify the reason why TeV gamma-rays 
have until now been detected only from SN1006. 



With the above motivation, we have observed 
RXJ1713.7-3946 with the CANGAROO imaging TeV 
gamma-ray telescope in 1998. Here we report the result 
of these observations. 



2. Instrument and Observation 

The CANGAROO 3.8m imaging TeV gamma-ray tele- 
scope is located near Woomera, South Australia 



(136°47'E, 31°06'S) (Hara et al. |1993[ ). A high resolution 
camera of 256 photomultiplier tubes (Hamamatsu R2248) 
is installed in the focal plane. The field of view of each 
tube is about 0°.12 x 0°.12, and the total field of view 
(FOV) of the camera is about 3°. The pointing accuracy 
of the telescope is ~ 0°.02, determined from a study of 
the trajectories of stars of magnitude 5 to 6 in the FOV. 
RX J1713.7— 3946 was observed in May, June and August 
in 1998. During on-source observations, the center of the 
FOV tracked the NW rim (right ascension 17 h ll m 56 s .7, 
declination -39°31'52".4 (J2000)), which is the bright- 
est point in the remnant in hard X-rays (Koyama et al. 
1997). An off-source region having the same declination 
as the on-source but a different right ascension was ob- 
served before or after the on-source observation for equal 
amounts of time each night under moonless and usually 
clear sky conditions. The total observation time was 66 
hours for on-source data and 64 hours for off-source data. 
After rejecting data affected by clouds, a total of 47.1305 
hours for on-source data and 45.8778 hours for off-source 
data remained for this analysis. 



3. Analysis and Result 

The standard method of image analysis was applied for 
these data which is based on the well-known parameter- 
ization of the elongated shape of the Cerenkov light im- 
ages using "width," "length" "concentration" (shape), "dis- 
tance" (locat ion), and the im age o rientation angle " alpha " 
(Hillas |l985| , Weekes et al. |l989[ Reynolds et al. |1993[ ). 
However, the emitting region of TeV gamma-rays in this 
target may be extended, as in the case of SN1006. For 
extended sources, use of the same criteria as for a point 
source in the shower image analysis is not necessarily op- 
timal. We made a careful Monte Carlo simulation for ex- 
tended sources of various extents and found the distribu- 
tion of the shower image parameter of width, length, and 
concentration for gamma-ray events is essentially the same 
within a statistical fluctuation as in the case of a point 
source. However, the simulation suggests that we should 
allow a wider range dependent on the extent of the source 
for the parameter of distance and alpha to avoid overcut- 
ting gamma-ray events. In this analysis, gamma-ray-like 
events were selected with the criteria of 0.°01 < width < 
0.°11, 0.°1 < length < 0.°45, 0.3 < concentration < 1.0 
and 0.°5 < distance < l.°2. 
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Figure [l]a shows the resultant alpha distribution when 
we analyzed the distribution centered at the tracking point 
(right ascension 17 h ll m 56 s .7, declination -39°31'52".4 
(J2000)), which is the brightest point in the remnant in 
hard X-rays (Koyama et al. 1997). The solid line and the 
dashed line indicate the on-source and off-source data re- 
spectively. Here we have normalized the off-source data to 
the on-source data to take into account the difference in 
observation time and the variation of trigger rates due to 
the difference in zenith angle between on- and off-source 
data and due to subtle changes in weather conditions. The 
value of the normalization factor (3 is estimated to be 1.03 
from the difference in total obsevation time for on- and off- 
source measurements. On the other hand, the actual value 
of the normalization factor (3 is estimated to be ~ 0.99 
from the ratio of N on /N Q s, where N on and N g indicate 
the total number of gamma-ray-like events with alpha be- 
tween 40° and 90° for on- and off-source data respectively. 
We selected the region with alpha > 40° to avoid any "con- 
tamination" by gamma-rays from the source, in the knowl- 
edge that the source may be extended. The small discrep- 
ancy in the two estimates of the value (3 might come from 
a slight change in the mirror reflectivity during the ob- 
servation due to dew. Here we adopt the value 0.99 for /3 
in the following analysis by taking the small discrepancy 
into the systematic errors due to the uncertainty in the 
mirror reflectivity as shown below. Figure 0b shows the 
alpha distribution of the excess events for the on-source 
over the off-source distribution shown in Figure 0a. A 
rather broad but significant peak can be seen at low al- 
pha, extending to 30°. The alpha distributions expected 
for a point source and several disk- like extended sources of 
uniform surface brightness with various radii centered at 
our FOV were calculated using the Monte Carlo method. 
These distributions are shown in the same figure. The 
alpha distribution of the observed excess events appears 
to favour a source radius of ~ 0°.4, which suggests the 
emitting region of TeV gamma-rays is extended around 
the NW rim of RXJ1713. 7-3946. The statistical signifi- 
cance of the excess is calculated by (iV on (a) — 0N o fi(a)) / 
\/N on (a) + f3 2 N s(a), where N on (a) and N g(a) are the 
numbers of gamma-ray-like events with alpha less than 
a in the on- and off-source data respectively. The signifi- 
cance at the peak of the X-ray maximum was 5.6a when 
we chose a value of alpha 30° considering the result of the 
Monte Carlo simulation shown in Figure 0b. 

In order to verify this extended nature, we examined 
the effects of the cut in shape parameters on the alpha dis- 
tribution by varying each cut parameter over wide ranges. 
We also produced alpha distributions for different energy 
ranges and data sub-sets. Similar broad peaks in the alpha 
distribution persisted through these examinations. Also 
we examined more recent data from PSR1706— 44 from 
July and August 1998 and obtained a narrow peak at al- 
pha < 15°, as expected for a point source. This confirms 
that the extended nature of the TeV gamma-ray emitting 




(b) 



20 40 60 80 
Alpha (degrees) 



20 40 60 
Alpha (deg 



80 

rees) 



Fig. 1. (a) Distributions of the orientation angle "alpha" 
for gamma-ray-like events with respect to the center of the 
field of view, which for on-source data corresponds to the 
NW rim of RX J1713. 7-3946. The solid line and dashed 
line are for on-source and off-source data respectively, (b) 
Distribution of the excess events of the on-source over the 
off-source level shown in Figure 0a, shown as the shaded 
bins. The vertical bars for several bins indicate plus and 
minus one standard deviation which is approximately the 
same for all bins. The expected alpha distribution for a 
point source (dotted line), and disk- like sources with a 
radius of 0°.2 (dashed line) and 0°.4 (thick solid line) cen- 
tered at the FOV by the Monte Carlo method are also 
shown. Here the curves are normalized to the actual ex- 
cess number of gamma-ray-like events with alpha < 30°. 



region does not come from some malfunction of our tele- 
scope system and/or systematic errors in our data analy- 
sis. A similar, but not as broad, alpha peak was seen for 
SN1006 (Tanimori et al. |l998b[) . 

In order to see the extent of the emitting region, 
we made a significance map of the excess events around 
the NW rim of RXJ1713. 7-3946. Significances for alpha 
< 30° were calculated at all grid points in 0°.l steps in 
the FOV. Figure shows the resultant significance map of 
the excess events around the NW rim of RX J1713.7— 3946 
plotted as a function of right ascension and declination, in 
which the contours of the hard X-ray flux (Tomida 1999) 
are overlaid as solid lines. The solid circle indicates the 
size of the point spread function (PSF) of our telescope 
which is estimated to have a standard deviation of ~ 0°.25 
for alpha < 30° based upon Monte Carlo simulations for 
a point source with a Gaussian function. The area which 
shows the highest significance in our TeV gamma-ray ob- 
servation coincides almost exactly with the brightest area 
in hard X-rays. The region which shows the emission of 
TeV gamma-rays with high significance (> 3<r level) ex- 
tends wider than our PSF and appears to coincide with 
the ridge of the NW rim that is bright in hard X-rays. It 
extends over a region with a radius of ~ 0°.4. This region 
persisted in similar maps calculated for several values of 
alpha narrower than 30° . 

The integral flux of TeV gamma-rays was calculated, 
assuming emission from a point source, to be (5.3 ± 0.9 
[statistical] ± 1.6 [systematic]) x 10~ 12 photons cm~ 2 
s" 1 (> 1.8 ± 0.9 TeV). The flux value and the sta- 



4 



H. Muraishi et al.: TeV gamma-ray emission from the SNR RX J1713. 7-3946 




Relative Right Ascension (degrees) 

Fig. 2. Contour map of significance around the NW rim of 
RX Jf 743.7— 3946 centered at the region brightest in hard 
X-ray emission (right ascension 47 h 44 m 56 s .7, declination 
-39°34'52".4 (J2000)). White lines indicate the signifi- 
cance level. The contours of the 0.5-40 keV band of the 
X-ray flux (from Tomida 4999) also are overlaid as solid 
lines. The solid circle indicates the size of our PSF. 



tistical error were estimated from the excess number of 
iVon(30°) - /3iV off (30°), where the value of 30° for alpha 
is chosen by the argument mentioned before. The causes 
of the systematic errors are categorized by uncertainties 
in (a) assumed differential spectral index, (b) the loss of 
gamma-ray events due to the parameter cuts, (c) the es- 
timate of core distance of showers by the Monte Carlo 
method, (d) the trigger condition, (e) the conversion fac- 
tor of the ADC counts to the number of photo-electrons, 
and (f) the reflectivity of the reflector. These errors from 
(a) to (f) are estimated as 45%, 22%, 3%, 12%, 40%, and 
8% for the integral flux and 24%, 2%, 8%, 20%, 29%, and 
17% for the threshold energy, respectively. The total sys- 
tematic errors shown above are obtained by adding those 
errors quadratically. 

To summarise, all our observed data support the hy- 
pothesis that the emitting region of the NW rim is ex- 
tended. In general, the value of the effective detection area 
of the telescope system for extended sources would be re- 
duced by some factor from that for a point source, because 
the gamma-ray detection efficiency decreases with the dis- 
tance of emitting points from the center of the FOV when 
we observe with a single dish. We calculated the efficiency 
as a function of the distance by the Monte Carlo method 
by analyzing the data with the same criteria as applied 
to the actual data. We estimated the value of the correc- 
tion factor to the effective area to be ~ 1.2 for our target 
by integrating the efficiency over the distance for an ex- 
tended disk-like source of uniform surface brightness with 
a radius of 0°.4. The factor of 1.2 is less significant than 
the systematic errors estimated above. 



4. Discussion 

The SNR RX J1713. 7-3946 is reminiscent of SN4006 both 
in the synchrotron X-ray emission from the shell far from 
the centre of the remnant and also in the TeV gamma-ray 
emission from an extended region coincident with that of 
the non-thermal X-rays. This suggests that the particles 
responsible for the emission of the high energy photons 
are accelerated in shocks. 

There are several possible emission processes of TeV 
gamma-rays: the emission induced by accelerated protons 
(by the ir° decay process) and by electrons - through 
bremsstrahlung and/or the Inverse Compton (IC) pro- 
cess. The expected integral flux of gamma-rays above our 
threshold energy of ~ 4.8 TeV by the 7r° decay process 
is estimated to be < 4 x 40~ 14 photo ns cm~ 2 s^ 1 (Drury 



et al. 4994 , Naito & Takahara 4994 ), where we assume 
the distance and the upper limit for the number den- 
sity in the ambient space of the remnant as 6 kpc and 
0.28 atoms/cm 3 , respectively (Slane et al. 1999). This flux 
value is too low to explain our observed flux, even taking 
into account the large uncertainties in the estimates of the 
distance and the ambient matter density of the remnant 
(Slane et al. 1999, Tomida 1999). However, there remains 
the possibility of some contribution of the 7r° decay pro- 
cess if the remnant is interacting with a molecular cloud 
located near the NW rim (Slane et al. 199E). The relative 
contribution in emissivity of the bremsstrahlung process 
compared to ir° decay process is estimated as ~ 10 %, 
assuming the flux ratio of electrons to protons is ~ 1/100 
and that bo th ha ve power law spectra with the index of 
2.4 (Gaisscr 1990 ), indicating this process is also unlikely 
to dominate. Therefore, the most likely process for TeV 
gamma-ray emission seems to be the IC process. 

Under this assumption, the magnetic field strength in 
the supernova remnant can be deduced from the relation 
L syn /Lic = Us/Uph between the IC luminosity Lie and 
synchrotron luminosity L syn , where Ub = B 2 /8ir and J7 p h 
are the energy densities of the magnetic field and the tar- 
get photon field, respectively. L syn and L\q in the above 
formula must be due to electrons in the same energy range. 
The value of L syn which should be compared with our TeV 
gamma-ray data is estimated from the ASCA result to be 
L syn = L ASCA f^E-^dE /f^E^dE, extrap- 
olating the synchrotron spectrum with the same power law 
out of the energy range of 0.5— 10 keV covered by ASCA 



(Tomida |l99S| ). Here Lasca = 2.0 x 40~ 10 erg cm" 2 s" 1 



is the X-ray luminosity in the 0.5—40 keV energy band 
observed by ASCA from the NW rim of the remnant and 
the power law index of —4.44 is the mean value for in- 



dex of X-rays in the same energy range (Tomida 1999). 
E™ 1 ^ = 0.14(f?/T0/zG) keV is a typical synchrotron pho- 
ton energy emitted by electrons which emit 1.8 TeV pho- 
tons (the threshold energy of our observation) by the IC 
process when we assume the target photons to be from the 
CMBR. The value of Lie is calculated to be 4.2 x lO" 11 
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erg cm -2 s -1 from our result for the number of photons 
of TeV gamma-rays, and using the fact that the spectra 
of synchrotron photons and IC photons follow the same 
power law when the electrons have a power law spectrum. 
Thus inserting L syn , Lie, an d ?7ph = 4.2 x 10 -13 erg cm -3 
of the energy density for the CMBR into the above re- 
lation, we can solve for the magnetic field strength B. 
Finally, the magnetic field at the NW rim is estimated to 
be ~ 10.8^G. The extrapolation used to estimate L syn is 
reasonable, because E™™ is estimated to be 0.15 keV; this 
is not so different from the minimum energy of the ASCA 
band (0.5 keV). 

The electrons responsible for the synchrotron and IC 
photon emissions are likely to have been accelerated by 
the shocks in the remnant as discussed above. If the max- 
imum electron energy is limited by synchrotron losses, this 
maximum energy can be estimated by equating the cool- 
ing time due to synchrotron losses with the time scale of 
acceleration by the first Fermi process in a strong shock 
as ~ 50(K/2000kms- 1 )(S/10AtG)-°- 5 Te V, wh ere V s is 
the shock velocity (Yoshida & Yanagita 1997). On the 



other hand, equating the acceleration time with the age 
of the remnant, the maximum energy can be expressed 
- 180(y s /2000km s- 1 ) 2 (B/10AtG)(t agc /10000year) TeV. 
In either case, whether it is synchrotron losses or the age 
of the remnant that limits the maximum electron energy 
(Reynolds & Keohane 1999), electrons should exist with 
energies high enough to emit the observed synchrotron 
X-rays and TeV gamma-rays by the IC process. 

It is notable that both RXJ1713. 7-3946 and SN1006 
have relatively low magnetic field strengths and low mat- 
ter densities in their ambient space. These common fea- 
tures may have arisen if the magnetic field was 'frozen 
in' to the matter without amplification other than by 
compression by shocks and may be the reason why elec- 
trons are accelerated to such high energies. These facts 
may also explain the radio quietness (Green 1998|) and 
the weak emissivity of ir° decay gamma-rays of the rem- 
nants. For SN1006, the low matter density in the ambient 
space might result from the remnant being located far off 
the galactic plane and the supernova being of type la. For 
RX J1713.7— 3946, the low matter density may be caused 
by material having been swept out by the stellar wind 
of the supernova progenitor (Slane et al. 1999). The low 
magnetic field and the low matter density in the ambient 
space of SN1006 and RXJ1713. 7-3946 may explain why 
TeV gamma-rays have been detected so far only for these 
two remnants. 

In conclusion, we have found evidence for TeV gamma- 
ray emission from RX J1713.7— 3946 at the level of 5.6 



sigma. If confirmed (a la Weekes 1999), this would be the 
second case after SN1006 to show directly that particles 
are accelerated up to energies of ~ 100 TeV in the shell 
type SNR. 
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